Thermodynamics of Ising Spins on the Star Lattice 
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There is a new class of two-dimensional magnetic materials polymeric iron (III) acetate fabricated 
recently in which Fe ions form a star lattice. We study the thermodynamics of Ising spins on the star 
lattice with exact analytic method and Monte Carlo simulations. Mapping the star lattice to the 
honeycomb lattice, we obtain the partition function for the system with asymmetric interactions. 
The free energy, internal energy, specific heat, entropy and susceptibility are presented, which can 
be used to determine the sign of the interactions in the real materials. Moreover, we find the rich 
phase diagrams of the system as a function of interactions, temperature and external magnetic 
field. For frustrated interactions without external field, the ground state is disordered (spin liquid) 
with residual entropy 1.522 . . . per unit cell. When a weak field is applied, the system enters a 
ferrimagnetic phase with residual entropy InA per unit cell. 

PACS numbers: 75.10.Hk, 75.30.Kz, 75.40.-s, 64.60.-i 



I. INTRODUCTION 

Spin systems with geometrical frustration have both 
fundamental and practical importance. Theoretically, 
lots of interesting phenomena have been found in the 
geometrically frustrated systems, like the antiferromag- 
netic triangular lattice, kagome lattice. The systems 
can remain disordered even at absolute zero tempera- 
ture because of the competitive magnetic interactions. 
For example, an antiferromagnetic triangular lattice has 
a residual entropy so = 0.3281... per unit cell 1 . The 
frustration effect has important application in achieving 
a lower temperatures through the adiabatic demagneti- 
zation compared with other methods. When the temper- 
ature, external magnetic field, and other factors are con- 
sidered, the geometrically frustrated systems can show 
very rich phase diagrams^ A typical case is that the mag- 
netocaloric effect can be enhanced near the phase transi- 
tion points when a finite external field is applied*^— Of 
considerable interest has been searching for geometrically 
frustrated systems. Some new frustrated materials have 
been fabricated and studied recently, such as Ro 2 Ti 2 7 , 
Dy 2 Ti 2 7 , and Cu 9 X 2 (cpa) 6 -xH 2 0(X = F, CI, Br)&~— 

In 2007, a new class of geometrically frustrated mag- 
netic materials polymeric iron (III) acetate — was fab- 
ricated, in which Fe ions form a two-dimensional lattice 
referred as star lattice. Experiment has found that the 
materials exhibit spin frustration and have two kinds of 
magnetic interactions: intratrimer Jt and intertrimer Jjj 
shown in Fig. [TJ The Fe ion has a large spin, which is 
S = 5/2. The system may have the paramagnetic ground 
state because of the geometrical and quantum fluctua- 
tions. 

The Ising model on the star lattice with uniform ferro- 
magnetic couplings has been solved^ The critical tem- 
perature K c = 0.81201 was given with K c = f3 c J. The 
Bose-Hubbard model on the star lattice has also been 
studied using the quantum Monte Carlo method.— Re- 
cently, the edge states and topological orders were found 
in the spin liquid phases of star latticed Even though 



the 5 = 1/2 quantum Heisenberg model is considered 
to be appropriate to help studying the new material be- 
cause of its quantum fluctuations in the ground state, 17 
the Ising model on the star lattice is still very important 
especially for the non-uniform case. In real materials, 
the Fe ion has 5 = 5/2 which is close to the classical 
limit and the magnetic system shows two types of inter- 
actions^ Therefore, it is important to study Ising spins 
on the star lattice with the asymmetric interactions, es- 
pecially for the frustrated case. 

In this paper, we aim at the thermodynamics of Ising 
spins on the star lattice with asymmetric interactions us- 
ing the exact analytic methods and Monte Carlo sim- 
ulations. We present the phase diagram as a function 
of interactions, temperature and external magnetic field. 
There is a clear difference between the Jd > case and 
the Jd < case. Our study provides useful information 
for determining the sign of Jd- 

This paper is organized as follows. The model is de- 
scribed in Sec. II. In Sec. Ill, we map the star lattice 
to honeycomb lattice and get the exact results. Sec. IV 
presents the phase diagrams as functions of interactions 
and external magnetic field. In Sec. V, the Monte Carlo 
results for the heat capacity and susceptibility are given. 
Finally, in Sec. VI we summarize the results. 



II. MODEL 

The structure of star lattice and its unit cell are illus- 
trated in Fig. [TJ We study Ising spins on the star lattice 
with two kinds of nearest-neighbor interactions, the in- 
tratrimer coupling Jt and intertrimer coupling Jd- The 
Hamiltonian is 



H = -J T ^ " (Jidj - J D ^ Gi><Jj> - A^a*, (1) 

<i'j'> i 



<ij> 



where < ij > runs over all the nearest neighbor spin 
pairs, Jt is the intra-triangular interaction and Jd is the 
inter-triangular interaction, h is the external magnetic 
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FIG. 1: (Color online). The star lattice. The dashed frame 
represents a unit cell of the star lattice. There are six spins 
per unit cell. 



field, and Oi = ±1. The unit cell of the star lattice 
contains six spins shown in Fig. [TJ If we use and 
Njj to denote the total numbers of Jt- and J^-bonds, 
we have Nt : Njj =2:1. The analytic result of partition 
function is obtained for h = 0. 

For simplicity, we use \Jt\ as the units of energy in 
the following. The corresponding phase diagrams are ac- 
tually in a three-dimensional parameter space, j^t^jj^j 

and tt-t. 

\Jt\ 



III. EXACT SOLUTION IN ZERO FIELD 

In this section, we study the exact analytic results of 
Ising model on the star lattice in zero magnetic field 
{h = 0). Using a sequence of A — Y transformation and 
series reductions^, we can transform the Ising model on 
the star lattice into one on the honeycomb lattice whose 
partition function has been exactly solved using the Pfaf- 
fian metho d 18 ! 19 . Besides the exact analytical results, we 
expand the partition function in series for some special 
cases. 
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FIG. 2: (Color online). Transformation of star lattice to hon- 
eycomb lattice, (a) Depicts a section of the star lattice with 
two couplings, Jd, Jt- By applying the A — Y transformation 
and then we obtain (b) . After that take the two bonds Jd , J\ 
in series and we obtain (c) where generates a new coupling Ji . 
Finally, take J\ and J2 in series and the honeycomb lattice 
(d) is obtained. 



change couplings of Fig. [2] can be written as 

1 



U 



t T + t' 1 - 1 

th = ht2- 

We write th in terms of tr and tn directly 

trio 



t h 



For convenience, we can rewrite it in terms of xi 



Xh 



xd + (2 + xd)x\ 
1 + (l + 2x D )x?r 



(2) 

(3) 
(4) 



(5) 



(6) 



A. Effective coupling on the equivalent honeycomb 
lattice 

The results of A — Y transformation and series re- 
duction are given in Ref. 10. Using the variables U = 
tanh/3Ji and Xi = e _2/3Ji , the relations among the ex- 



B. Phase boundary 

It is known that the critical temperature of the hon- 
eycomb lattice Ising model is given by x c h = 2 - a/3. 20 
Having mapped star lattice to honeycomb lattice, we can 
substitute this into the equivalent coupling in Eq. (j6j). 
Thus, an implicit equation for the critical temperature 
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FIG. 3: (Color online). Phase diagram of the star lattice Ising 
model in the (Jd,T) plane, for J T = 1 and h = 0. The thick 
curve is the exact solution. It illustrates that when Jd > 0, 
the phase is ordered. The ordered phase is ferromagnetic. On 
the contrary, when Jd < the phase is immediately disor- 
dered (paramagnetic). When Jt — — 1, the phase diagram 
is below T-axis (not shown in the figure), which implies that 
there is no phase transition in this phase. 



of star lattice Ising model can be obtained as, 



-2(3 C J D 



(2- 



2{3 c J D y-4{3 c J T 



= 2-^3 (7) 



This result is plotted in Fig. 03 When Jd is ferromagnetic 
(Jd > 0) and strong enough, the critical temperature 
saturates at a finite value, i.e. T c /\Jt\ ~ 4/Zn[3/(2\/3 — 
3)] w 2.14332. The critical temperature drops to zero 
as Jd approaches zero. When Jt ~ \Jd\, the curve is 
approximately linear with T c /\Jt\ ~ 1.23151. 

Furthermore, when xd = %t = x, Eq.([6]) reduces to 
the result of star lattice with equivalent couplings^ We 
get x c = 0.19710, or equivalently, K c = 0.81201. 

Since all the factors obtained here are analytical, the 
singularity in the partition function remains when we 
transform the star lattice into the honeycomb lattice. 
The phase transition is the same as the honeycomb lat- 
tice where a continuous second-order transition happens. 

If Jd is antiferromagnetic (Jd < 0), we get the neg- 
ative critical temperature, which implies no phase tran- 
sition existing in this case. It can be used as a criterion 
for experimentalists to determine if the couplings in a 
real material is ferromagnetic or antiferromagnetic. If 
one finds a phase transition in the real material, we can 
conclude that the couplings Jd and Jt should be both 
ferromagnetic. In there is no long range order found, it 
means that at least one kind of nearest neighbor cou- 
plings is antiferromagnetic in the material. 



C. Partition function 

Since we have utilized the A — Y transformation and 
series reductions to map the star lattice to a honeycomb 
lattice, the partition function per unit cell, z 3 , of the 



star lattice Ising model is equivalent to that of the hon- 
eycomb lattice zh multiplied by some coefficients which 
result from the transformation. These coefficients are as 
follows, 




z 2 = (1 + XiX D ) 



Z 3 = (1 + XiX 2 ) 



X\X2 



(8) 
(9) 
(10) 



Therefore, the total partition function of the star lattice 
is 



(ii) 



where zh is calculated using the Pfaffian method^ We 
rewrite it here, 



= — exp{±n[w(x h )]} (12) 



where 

n(w) = 

and 



r 2n dp r 2 

Jo ^Jo 



dq 
2^ 



\n(w — cosp — cos q — cos(p + q)) 

(13) 



1 - 2x h + 6x h 2 - 2x h 3 + x h A 

w{Xh) = 2x h{ i-x h r • (14) 

We can rewrite Q(w) and get a more accurate numerical 
evaluation according to the singularities of the integrand. 



2 r /2 

Q(w) = — I dp In 



cosp + arccosh 



w — cos 2p 



2 cosp 



(15) 



The partition function of the star lattice Ising model is 
therefore 

z s (x T ,x D ) = ^(x T ,x D )ex.p [\Vt(w(x h (x T ,x D )))\ (16) 
where 

S&(x T ,x D ) = x^ 3 x D ^(l - x T )(l - x 2 D ) 

x W2(l + x\ + 2x D x T )(x D + 2x\ + x D x T ) (17) 



The total partition function is given by Z s = , where 
N is the spin number of unit cell. Since the partition 
function is obtained, the internal energy, specific heat, 
entropy and free energy can be calculated from it. 



D. Energy 

Taking derivation of the partition function, the energy 
per unit cell of the star lattice Ising model can be ob- 
tained. 



u = 



d\nz 

d(3 = 

=D,T 



dxT dhiz dxD d In z 
df3 dxr df3 dxo 
<91n\I/ dxh dw dVt 

2 ^ 

OXo 



dxi dxh dw 



(18) 



where ^ is expressed in terms of the complete elliptic 
integral of the first kind, K^l 

dVt o 
dw 
xK 



7T(-W-l) 3 / 4 (-W + 3) 1 /4 

( 1 | w 2 -3 
V 2 



2(>+l)(-w-l)i/ 2 (-w+3) 



m) ■ (19) 



The plots of energy in units of \Jd\-, jj^j is illustrated 
in Figs H] and [5] for the unfrustrated case and frustrated 
case respectively 



E. Specific heat 




FIG. 5: (Color online). Thermodynamic functions for one 
unit cell vs temperature T for frustrated coupling Jd — —0.5 
and Jt = —1- The specific heat (red line) is no longer di- 
verging. Entropy (orange line) remains a none-zero value at 
T = and 6/n2 at high temperature. 



However, it arises again when R > 6 In the unfrustrated 
case, the sharp peak vanishes which implies no phase 
transition, consistent with the conclusion drawn from the 
phase diagram. 



By further derivation, c = J^, the heat capacity per 
unit cell can be obtained. The details are shown in 
Ref. |To|- Here we show the plots of c in Figs |4] and [5] 

In the unfrustrated case, the specific heat c shows a 
sharp peak at T w 0.74|Jx| where a phase transition 
happens. The phase transition point is consistent with 
the result of Eq. (J6j). In addition, there is a broad hump 
at higher temperature because of the flopping of spins. 
Moreover, this hump changes with R = \ jP-\. It is ob- 
vious when R < 1 and becomes indistinct when R = 1. 




T 
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FIG. 4: (Color online). Thermodynamic functions for one 
unit cell vs temperature T for the unfrustrated case Jd — 0.5 
and Jt = 1. The specific heat (red line) diverges as T % 
0.74|Jt| revealing that there is a phase transition from the 
ferromagnetic phase to the parametric phase. Energy (blue 
line) is shown as —u(T). Entropy (orange line) approaches 
zero when T — ► and saturates at 6/n2 as T — ► oo. 



F. Zero-temperature limit: residual entropy 

The plots of entropy are shown in Figs. |4] and [U 
Nonetheless, we can expand the partition function in se- 
ries to gain more information about the residual entropy 
in the low temperature limit. 

In the case of Jd > 0, the partition function can be 
expanded as 

3 3 
Inz = —-InxD — SIuxt + -x D + ... (20) 

u = -3\J D \ - 6\J T \ + 6| J D \e~ m Jd1 + ... (21) 

The residual entropy is therefore when T — 0. 
However, when Jjj < 0, the model becomes frustrated. 
In this way, when T — >• 0, j3 — >• oo, which means 
xd,xt oo, w — >• oo. Therefore, lnQ(w) becomes 
~ ln(w). Expanding ln(z), we get 

1 3 

Inz = -ln21 + -Iuxd + Inxr + ... (22) 

u=(3\J D + 2\J T \) + ... (23) 
These results contribute to the residual entropy by 

s lim lim (In Z + f3u) = -ln21. (24) 

(3 J — y — oo f3 J d — y — oo 2 

Thus, the frustration of the system leads to a ^ln21 w 
1.522 residual entropy per unit cell when T 0. One 
can confirm that, this value is consistent with the en- 
tropy at T = in Fig. [5] The residual entropy per site 
is approximately 0.254, smaller than the triangle lattice, 
TKL, and kagome lattice ^ 1Q i 22 . Therefore, the star lat- 
tice is less frustrated compared to them. 
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IV. PHASE DIAGRAMS AT ZERO 
TEMPERATURE 

In this section, we present the phase diagrams at zero 
temperature along with some thermodynamic properties 
such as energy, magnetization and entropy. By calculat- 
ing the ground state energy of the star lattice, we derive 
the full phase diagram for the system. Since the phase di- 
agram at zero field is already shown in Fig. EJ we focus on 
the none-zero field case in this section. The correspond- 
ing results are summarized in Figs. [TUI and fTTI according 
to the sign of Jt- 



A. Zero Field (Phase V and VI) 

The phase diagram for zero field as a function of cou- 
plings is showed in Fig. 02 When Jd > 0, the phase is 
ordered and ferromagnetic. When Jd < 0, the phase is 
frustrated with a residual entropy sq = |/n21. When 
Jt = —1, the phase is located in the negative section of 
T-axis , which reveals that there is no phase transition 
in this phase. The disordered and ordered phases are 
labeled by V and VI in Figs. [TOland fTTI respectively. 

In phase V, the system is fully frustrated. We find 
18 degenerate ground states for each unit cell. However, 
as shown in Eq. (f2l|) , the residual entropy is not InlS 
but \ln2\. This is similar as the triangular lattice whose 
residual entropy can not be obtained by counting the 
number of ground states in a unit cell. U 




FIG. 6: (Color online). Spin configurations of the degenerate 
states of phase II. The two spins are coupled by Jt- Each 
triangular has exactly two spins pointing up. 




FIG. 7: (Color online). Spin configurations of the degenerate 
states of phase III. The two spins are coupled by Jt- Each 
triangle has exactly two spins pointing up. 




FIG. 8: (Color online). Spin configurations of the degenerate 
states of phase IV. All the configurations have only one spin 
points down. 



B. Saturated ferromagnetic phase (Phase I) 

When the external field is strong enough, i.e. h > 
Max{2\ Jt|, I Jd\ + 2|Jr|}, the phase is ferromagnetic 
where all spins are lined up. It is obvious that this state 
has energy u = —Jd — 6Jt — 6/i, magnetization m = 6 
and entropy s = per unit cell. 



C. Phase II 

When the field is weaker, e.g. < h < 2\Jt\ and Jt < 
0, Jd > 0, it is a phase with u = — Jd + 2 Jt — 2h, m = 
2, s = Inb. The spin configurations of the degenerate 
ground states of this phase are shown in Fig. [6j The two 
spins connected by Jd become aligned due to the positive 
Jd and the weak field h. 
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FIG. 9: (Color online). Spin configurations of the degenerate 
states of phase VII. This is a interesting phase because spins 
in the same triangle point the same direction. It can map to 
a honeycomb lattice with equivalent antiferromagtic coupling 
with high spins. 



D. Phase III 

If < ft < 2, the system is in a frustrated phase. 
We find four degenerate ground states in this phase con- 
tributing to the residual entropy s = InA. The other 
properties are given by u = Jd + 2 Jt — 2ft, and m = 2. 
The spin configurations are shown in Fig. In this case, 
the two spins connected by Jt become antiparallel since 
Jd is antiferromagnetic. 



E. Phase IV 

In the case of 2 < ft < |Jt>|+2|J t | and J T > 0, Jt> < 0, 
phase III evolves into phase IV, which has m = 4 and 



\ _h 

\ \Jt\ " 






m = 4 \ 

iv s=ln2 \ 


I 


m = 6 

5 = 


III m = 2 


II 


m — 2 
s = ln5 









-2 



\Jt\ 

FIG. 10: (Color online). Phase diagram of the star lattice 
Ising model in the (Jd, ft) plane for Jt > and T = 0. The 
phase diagram is symmetric under a sign change of ft. 



s = ln2. Only one spin points down in this phase and 
it should be one of the two connected by Jd- The spin 
configurations are shown in Fig. [U 

F. Phase VII 

Phase VII is a new phase when Jt becomes positive 
in none-zero field. In this phase, ft < \\Jd\ and Jd < 
0, Jt > 0, the spins on the same triangle are parallel, 
however, antiparallel to the neighboring triangles for the 
positive Jt and negative Jd- If we treat the three spins 
on the same triangle as a higher spin located in the center 
of the triangle, it becomes an antiferromagnetic phase in 
a honeycomb lattice. This state gives u = Jd — 6 Jt, s = 
ln2 and m=0. 



G. Phase diagram 

According to the discussion above, we now combine 
all the results together to obtain a full phase diagram. 
Fig. [10] show the phase diagram when Jt is antiferro- 
magnetic and thus Jt < and Fig. [TTJ on the contrary, 
shows the case when Jt is ferromagnetic. The phase di- 
agram is symmetric under the sign change of ft. 



V. MONTE CARLO SIMULATIONS 

In this section we show the Monte Carlo (MC) simula- 
tion results of the star lattice Ising model with different 





IV \ \J T \- 
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FIG. 11: (Color online). Phase diagram of the star lattice 
Ising model in the (Jd, ft) plane for Jt < and T = 0. 
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combinations of parameters, which helps to corroborate 
our analytic predictions. Meanwhile, they allow us to 
calculate the magnetization and susceptibility at finite 
temperature. 

We choose system size L = 8, 16, 32, where L is the 
length of the unit cell for the star lattice, which means 
that the total number of spins N is N = 6L 2 , as there 
are six spins in each unit cell. The periodic boundary 
condition is used for the simulations. 

The specific heat c and magnetic susceptibility x dur- 
ing the MC simulations can be calculated using the 
fluctuation-dissipation theorem 



X 



(H 2 ) - (Hf 
NT 2 

(m 2 ) - (My 

7VT 



(25) 
(26) 



where (H) and (M) are respectively the Monte Carlo 
averages of the total energy (i.e., the Hamiltonian) and 
magnetization. 

Fig. [T2l shows the temperature dependence of heat ca- 
pacity per site at h — for a typical unfrustrated case 
Jd/Jt — 5, Jt > 0. The MC results are consistent with 
the exact analytic results. 

We also calculate the susceptibility from the MC sim- 
ulations, which is shown in Fig. [131 As L increases, the 
peaks become sharper and sharper, indicating a phase 
transition. 

We also study two different combinations of interac- 
tions for the unfrustrated case, which is Jt > 0, in 
Figs. [13] and [TU The temperature dependence of sus- 
ceptibility is found to be sensitive to the sign of Jd- 
When Jd > 0, the susceptibility has a sharp peak at 
the critical point between the ferromagnetic phase and 
paramagnetic phase. However, when Jd < 0, there is 
no such peak, which implies no phase transition in this 




FIG. 13: (Color online). Temperature dependence of suscep- 
tibility from the Monte Carlo simulations with L = 8, 16 and 
32 for J D = 5J T = 1. 





FIG. 14: (Color online). Temperature dependence of suscepti- 
bility from the Monte Carlo simulations with L — 8, 16 and 32 
for Jd = 5Jt = — 1. There is no apparent peak found, which 
means there is no continuous phase transition for Jt < 0. 
The exact solution in Fig. [3] gives the same result in this case. 



case, just as what we get in exact solutions. The shape of 
the susceptibility peak depends on the size of the system 
when Jd > 0, whereas for Jd < 0, the size of the system 
has no influence on the susceptibility. 



2 4 6 8 10 

FIG. 12: (Color online). Temperature dependence of heat 
capacity per site from exact solution (red line) and Monte 
Carlo simulations with L = 8, 16 and 32 for Jd — 5Jt — 1. 
The critical temperature T c « 2.1|Jt|, consistent with Fig. [3] 



VI. CONCLUSIONS 

In summary, we have studied the Ising model on the 
star lattice with two different exchange couplings Jt and 
Jd using both analytical method and Monte Carlo sim- 
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ulations. We have presented its thermodynamic proper- 
ties including internal energy, free energy, specific heat, 
entropy and susceptibility in the zero field. The phase 
transition temperature for Jt = Jd is exactly same as 
the one found in Ref^. There is no phase transition 
found if one of the couplings is antiferromagnetic. More- 
over, we have obtained the rich phase diagrams in terms 
of Jt, Jd and h at zero temperature. Monte Carlo sim- 
ulation is used to confirm the exact results and calculate 
the susceptibility. 

In the fully frustrated case, the residual entropy of the 
system can be expressed as a closed form (s = \ln2\ as 
showed in Eq. (f24|) which is consistent with the triangu- 
lar and Kagome lattices. The system is less frustrated 
compared to the other triangulated lattices. 

Our study provides a benchmark calculation for the 
thermodynamics of Ising spins on the star lattice, which 



can help experimentalists to investigate the real materi- 
als. 
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